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We present a simple paradigm for detection of an immobile target by a space-time coupled random 
walker with a finite lifetime. The motion of the walker is characterized by linear displacements at a 
fixed speed and exponentially distributed duration, interrupted by random changes in the direction 
of motion and resumption of motion in the new direction with the same speed. We call these 
walkers ’’mortal creepers”. A mortal creeper may die at any time during its motion according to an 
exponential decay law characterized by a finite mean death rate u) m . While still alive, the creeper 
has a finite mean frequency uj of change of the direction of motion. In particular, we consider 
the efficiency of the target search process, characterized by the probability that the creeper will 
eventually detect the target. Analytic results confirmed by numerical results show that there is 
an cu m -dependent optimal frequency w = oj op t that maximizes the probability of eventual target 
detection. We work primarily in one-dimensional (d = 1) domains and examine the role of initial 
conditions and of finite domain sizes. Numerical results in d = 2 domains confirm the existence of 
an optimal frequency of change of direction, thereby suggesting that the observed effects are robust 
to changes in dimensionality. In the d = 1 case, explicit expressions for the probability of target 
detection in the long time limit are given. In the case of an infinite domain, we compute the detection 
probability for arbitrary times and study its early- and late-time behavior. We further consider the 
survival probability of the target in the presence of many independent creepers beginning their 
motion at the same location and at the same time. We also consider a version of the standard 
“target problem” in which many creepers start at random locations at the same time. 

PACS numbers: 02.50.Ey, 05.40.Fb, 82.20.Pm 


I. INTRODUCTION 

There are many processes in nature in which one or more randomly moving entities detect one or more targets. 
Among many examples we mention binding or trapping processes in molecular environments, diffusion-limited chemical 
reactions, and the search for nutrients in a biological context such as predator-prey situations in ecological systems. 
In many of these “target problems” the question of interest is the efficiency with which the mobile entity finds the 
target as measured by the time it takes to accomplish this outcome. The mathematical quantities to be calculated for 
this purpose are first passage time distributions or their moments, especially the mean first passage time (MFPT) for 
an entity to reach the target. The above problem is then tantamount to the calculation of the survival probability of 
the target. Many such problems for moving entities that execute random walks or, essentially equivalently, diffusive 
motion, have been presented in ref. fl[ • Recently, the role of parallel and intermittent search strategies in the above 
context has also attracted considerable interest (see e.g. refs. i ID- 

Other types of motion have frequently been considered in the context of target capture, especially in the ecology 
literature [4, including Levy flights and Levy walks (see e.g. chapters 7 and 8 in ref. [5j and a recent review on 
Levy walks Q). Levy flights are characterized by instantaneous steps whose lengths are chosen from a probability 
distribution that decays as a power law. This leads to scale-free superdiffusive motion. Levy walks differ from Levy 
flights in that the steps are not instantaneous but occur at a constant speed. Other models of superdiffusion are based 
on step-step correlations 3 and have also been considered in the target capture problem (see e.g. the example with 
fractional Brownian motion cited in Q). Whereas normal diffusion is associated with a mean square displacement 
that grows linearly with time, ( x 2 ) ~ t, in superdiffusive motion the growth is superlinear. At the other extreme lie 
subdiffusive models, ones in which the growth is sublinear. These are often associated with crowded environments 
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where motion is impeded by traps, barriers, and dead ends. Recent work has extended target capture models to 
subdiffusive motions arising from continuous time random walks (CTRWs) [Ml- 

Stepping back from the outcomes of various groups of models in terms of the mean square displacement to the 
mesoscopic basis for these outcomes, we can broadly characterize the random motions by two distributions: the 
distribution of stepping times (in the language of continuous time) and the distribution of stepping lengths (in the 
language of continuous space). Normal diffusion appears if the first moment of the waiting time distribution and the 
two first moments of the stepping length distribution are finite. The ubiquitous example is an exponential distribution 
of stepping times and a Gaussian distribution of stepping lengths. Subdiffusive motion can be achieved with stepping 
time distributions that lack a first moment, that is, there is a strong probability of a long sojourn at one location 
before moving on, so much so that the mean stepping time diverges. Superdiffusive motion, on the other hand, is 
associated with distributions of stepping distances that lack a second moment, that is, the variance of the step length 
diverges. We note that in recent years work on target capture has been extended to more complex processes such as, 
for instance, ones displaying non-trivial spatial or space-time correlations. One such example is given by the capture 
of an evasive prey which is able to detect predators within a certain range and step back accordingly [131 ]. 

A feature of our own recent work is the recognition that in many problems the moving entity has a finite lifetime 
dill. In general, this finite lifetime is randomly distributed. We have considered a variety of lifetime distributions, 
from exponential to algebraic. The former is the distribution associated with a “unimolecular” concentration decay 
process entity —> 0, while the latter may be indicative of more complex, encounter-controlled kinetics leading to death 
or inactivation [22j]. 

Returning to the target problem, the goal of a great deal of work, whichever the underlying model of motion and 
death, is to establish if optimization is possible. The goal may be to maximize the probability of encountering the 
target (lowest survival probability) or to minimize it (greatest survival probability). The latter case is for instance 
relevant in ecology, e.g., if a moving prey is trying to avoid a predator. Apart from the aforementioned details of the 
stepping time and the stepping length distributions, the relevant parameters with respect to which the search process 
must be optimized vary with the details of the model, for example initial conditions, boundary conditions, size, or 
dimensionality of the system. In particular, the literature refers to the target problem when there are many walkers 
starting their walk at random locations. There is a vast literature on the target problem [3, HJ, but not on the target 
problem with mortal rather than immortal walkers. However, in recent days the literature on the latter topic and 
related problems is increasing rapidly fl5l - |2ll |H, |24| . 

In any case, interesting effects arise from the interplay of the lifetime of an entity and the time it takes it to reach 
the target. If the lifetime is much shorter than the time it would take an immortal entity to reach the target, then the 
problem is in a sense uninteresting because the target is simply hardly ever reached. In the reverse case, the fact that 
the entity has a finite lifetime does not matter because it most frequently reaches the target before it dies. Clearly, 
the most interesting scenarios arise when these times are comparable. Since these times are both distributed, the 
most interesting cases are those where there is a significant overlap of the two distributions. 

In this paper we examine the search efficiency of mortal Levy walkers whose step duration and step length are 
exponentially distributed. More explicitly, the motion of our walkers is characterized by ballistic displacements at a 
fixed speed , interrupted at random times by random changes in the direction of motion and resumption of motion in 
the new direction with the same speed as before. Following Hughes @1 1 we call these walkers “creepers” to distinguish 
them from “leapers”, that is, walkers performing instantaneous jumps. In addition, we assume that our creepers are 
mortal, i.e., they have a finite lifetime which we assume to be exponentially distributed here. 

In investigating the survival probability of the target and the optimization (in this work, minimization) thereof, 
the interplay of two parameters is of special interest, namely, the average death rate w m and the average frequency of 
reorientation w of the creeper. We also emphasize an important point of the model: the creeper may reach the target 
at any time during its trajectory, which need not coincide with the end of a linear displacement. Most (but not all) 
of our paper deals with this problem in one dimension, where reorientation is simply a direction reversal. This case is 
admittedly considerably simpler than the problem in higher dimensions, but we do show numerically that the results 
are qualitatively similar in two dimensions. 

We organize our presentation as follows. In Sec. El we set forth the formalism for the motion of a single mortal 
creeper as well as for a collection of mortal creepers. In Sec. EH we use this formalism to calculate the survival 
probability of a target in the one-dimensional case, both in an infinite and in a finite domain. In the former case we 
give an explicit expression valid for arbitrary times. Section HVl presents a brief numerical study of the two-dimensional 
problem. Finally, in Sec. [V]we summarize our main findings and possible extensions of this work. 
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II. GENERAL FRAMEWORK 

Consider a system consisting of an immobile, impenetrable hyperspherical target of radius R and a randomly 
moving creeper of negligible spatial extent initially located at a certain distance from the target surface. The target 
is detected by the searcher as soon as the latter intersects the target surface. Our goal is to calculate the survival 
probability of the target. 


A. Master equation describing the motion of the creeper 

The motion of the creeper in one dimension starting from initial position Xq is described by the equation 

nt roo 

j(x,t\x o) = / dt' dx'jix-x^t-t'lxo^ix^t^+jo. (1) 

Here j(x, t\xo)dx dt is the probability that the creeper completes a displacement and thus performs a turn (direction 
reversal) during the infinitesimal time interval [t,t + dt] while dwelling in the infinitesimal segment [x,x + dx], and 
tR(rc, t) is the probability density of having performed a displacement a; in a time interval t between consecutive turns. 
Thus, Eq. (JT]) represents a detailed bookkeeping of the possible ways in which a creeper may finish a displacement at 
a given time in terms of all the possible positions and times of the previous turn. The deterministic initial condition is 
embodied in the choice jo(x) = 8(x — xo)8(t). In quantities that are conditional on the initial condition the dependence 
on initial time t = 0 is understood, e.g., j(x,t\xo) is an abbreviated notation for j(x, t\xo, 0). 

The description of the motion given by m is complemented by the equation 

nt /-OO 

p(x,t\xo) = / dt' / dx'j(x — x', t — t'\xo)4>(x', t r ), (2) 

J 0 J —oo 

which relates the statistical properties of the location of the creeper at time t to those of its sojourn in different regions 
of space for different periods of time. The quantity p(x, t\xo) dx is the probability of finding the walker inside the 
interval [x, x + da;] at time t provided it was at Xq at time t = 0, and 4>{x, t ) dx is the probability that a step remains 
unfinished after a time interval t since its start, during which time a displacement x was performed (if the creeper 
velocity is Vf > 0, then obviously oc 8(x ± v/t) - see below). 

The particular case of a walk with an exponentially decaying waiting time pdf is described by the following joint 
density: 


i&(x,t) =— [6(x + Vft) + d(x — Vft)]u>e ut , (3) 

where, as introduced earlier, w is the average rate at which displacements are terminated by a turn, implying that 
w -1 is the mean duration of a displacement in one direction. We also call this the characteristic “persistence time” 
of the creeper. The product form of ^(a^t) reflects the fact that the speed of each displacement is independent of its 
duration. For simplicity, we assume that the beginning of the first displacement coincides with the initial time t = 0, 
and that the creeper initially has equal probability to move to the left or to the right. The (5-functions in fact indicate 
the constant modulus Vf of the creeper’s velocity. Given the density in Eq. J3]), the corresponding expression for the 
density <j>(x, t) is obtained by computing the probability f f °° lo e~ ut dt' that a given step is unfinished after a time t: 

<j)(x, t) = i [8(x + Vft) + S(x — Vft)\ I to e~ ut dt' = i [<5(a; + Vft) + S(x — Vft)\ e~ ut = ^ X ’ ^ . (4) 

2 Jt 2 w 

Equations m provide a full description of the motion of an immortal creeper. 

We now introduce a finite lifetime for the creeper. We take the death process to be independent of the transport 
properties. The sojourn pdf p*(x,t\xo) for a mortal creeper can then be expressed in terms of the pdf p(x,t\xo) for 
an immortal creeper as 


p*{x,t\x 0 ) = ip{t)p(x,t\x 0 ), (5) 

where the mortality function <p(t) is the probability that the creeper has survived for a time interval t since it started 
moving, and </?(0) = 1. Throughout the remainder of this work we will use asterisks to denote quantities describing 
mortal creepers. 
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B. Probability of target detection 

We next turn to the calculation of the survival probability of a target centered at Xt g , that is, the probability 
that it has not been detected by a creeper. In the one-dimensional case we may assume without loss of meaningful 
information that the target radius is vanishingly small, R = 0. We first write the relevant equations for the case of 
immortal creepers and then show how to extend them to the case of mortal creepers. 


1. Single immortal creeper 


We introduce q(xt gi t\x', t') as the probability per unit time for the creeper to detect the target (that is, the target 
detection rate) at time t no matter how often this detection has occurred before, given that the creeper is found at x' 
at time t' < t. In particular, q(x tg ,t\xo, 0) = q(x tg ,t |xo) is the detection rate given that the creeper started its search 
at xo at the initial time t = 0. To detect the target, the creeper must step on location Xt g . Therefore q(xt g ,t\xo) is 
clearly proportional to the probability density p(xt g ,t\xo) introduced in the previous section (the special case with 
Xo = x tg will be dealt with separately). In what follows we assume perfect target detection, that is, the creeper 
detects the target with unit probability whenever it steps on location Xt g . 

We are ultimately only interested in the first encounter of the creeper with the target. The rate f(xt g ,t |xo) of 
first detection (probability per unit time to detect the target for the first time) for a creeper that starts at xo can be 
computed as a function of the overall detection rate q(xt g ,t\xt g ,t') from the renewal equation [26}: 


q(x tg ,t\x 0 ) = f(x tg ,t\x 0 ) + / q{x tg ,t\x tg ,t')f{x tg ,t'\x 0 )dt'. (6) 

Jo 

Here the rate at which the target is detected at time t has been decomposed into two mutually exclusive contributions. 
One contribution arises from first detection events at time t 1 given by the first term on the right hand side. The second 
contribution on the right hand side accounts for revisitations to the target location, that is, first visitations of the 
target at time t’ < t followed by a return to the target at time t (following any number of additional revisitations in 
between). Hence the first detection rate also appears in this term. 

Since the system is invariant with respect to a time shift, one has that p(xt g ,t\x' ,t') = p(xt g ,t — t'\x'), and thus 
q{xt g ,t\x', t') = q(xt g ,t — t'\x'). One therefore obtains the closed equation 

q(x t g, t\xo) = f{x tg ,t\x o) + / q(x t g,t - t'\x t g)f{x t g, t'\x 0 ) dt '. (7) 

Jo 

This is especially convenient because the integral is now a convolution so that the equation is greatly simplified in 
Laplace space. The Laplace transform of Eq. 0 readily leads to 


f(x tg ,s\x 0 ) 


q{x tg ,s\x 0 ) 

1 + q(xt g , s\x tg ) 


( 8 ) 


Note that we use the same notation for any function of time and its Laplace transform. The choice is made clear by 
the argument. The mean first detection time of the target (hereafter called the mean first passage time, MFPT) can 
then be obtained from this result for the Laplace transform of the probability per unit time of first detection, cf. [13 ■ 
We focus instead on the closely related survival probability of the target, which we assume to be killed as soon as 
the creeper arrives at the target location for the first time. The survival probability of the target is then given by 


In Laplace space, one then has 


S{t) = 1 - [ f{x tg ,t'\xo)dt'. 
Jo 


S{s) = - - - f{x tg ,s\x 0 ). 
s s 


In particular, the asymptotic survival probability Soc = S(t —► oo) is then given by 


See = 1 - 


POO 

/ f(x tg ,t\x 0 )dt = 1 - lim f(x tg ,s |x 0 ). 

Jo s “>° 


( 9 ) 


( 10 ) 


( 11 ) 


Since the searcher is immortal and the embedding geometry is one-dimensional, the motion is recurrent, that is, every 
location is guaranteed to be visited. This ensures that the target will eventually be detected with certainty, i.e., the 
integral on the right hand side of Eq. m is unity. These arguments can easily be generalized to the case of an 
arbitrary initial probability distribution p(x o) by taking the corresponding average in Eq. (1111) . 
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2. Single mortal creeper 


Equation (0 must be modified for mortal creepers since detection is now conditional on the survival function ip(t) 
of the creeper. To take this into account, we multiply Eq. © by tp(t) and rewrite the resulting equation as follows: 

q*(x tg ,t\x 0 ) = f*(x tg ,t\x 0 ) + [ q*(xtg,t\x tg ,t')f*(xtg,t'\x 0 )dt'. (12) 

Jo 

Here q*{x t g,t\x t g,t') = [ip(t)/(p(t')]q(xt g ,t\x t g,t') and f*(x tg ,t\x 0 ) = [ip{t)/ip(0)] f(x t g,t\x 0 ) = ip(t) f(x tg ,t\x 0 ). The 
quantities q* and /* have the same interpretation as q and /, but now weighted with the appropriate conditional 
survival probabilities. In particular, the ratio tp(t)/ip(t') is the survival probability of the creeper in the time interval 
[t',t\ conditional on having survived up to time t'. 

Since ip(i)/ip(t') depends on both t and t! separately and in general not on the time difference t — t\ the death 
process in general destroys the invariance of the propagator with respect to a time shift. It is then not possible to 
obtain a closed equation of the form ([7]) for the quantities with an asterisk. However, if the survival function decays 
exponentially, (p(t) = e~ Umt , one has ip(t)/ip(t') = ) anc j ^he shift invariance is preserved E2- For this 

particular case, one does find a renewal equation similar to ©, namely, 

q*{x tg ,t\x 0 ) = f*(x tg ,t\xo) + [ q*(xtg : t-t'\xtg)f*(xtg,t , \xo)dt'. (13) 

Jo 

In analogy to Eq. ([8jl. we find in Laplace space f*(x,s\xo) = q*(xt g , s|a;o)[l + q*{xt g , s\xt g )}~ 1 , and applying the shift 
theorem for the Laplace transform we obtain 


/•*/ I \ _ q{xtgi S -\- iX m |xo) _ .. | \ 

/ (^tgjSpo) — ill' I I -1 — f (Xtg, s ~h UJ m \x o). 

l + q(x t g,s + uj m \x t g) 

The survival probability at time t in the presence of an exponential mortality function then is 


s*(t) = l- f r{x tg ,t'\ Xo )dt', 

Jo 


whose Laplace transform yields 


S*(s) = 


1 - f*{Xtg,s\x 0 ) 


The asymptotic value of the survival probability is formally similar to that found in the absence of mortality, 


^ = 1 - 


/•OO 

/ f*{x t g,t\x 0 )dt = 1- limf*(x t g,s\x 0 ), 

Jo s-> - 0 


but it no longer vanishes. In particular, we now obtain the final result 

s * =l _ g {x tg ,s = w m \x 0 ) 

1 + q{x t g,S = Wm\x t g) 


(14) 

(15) 

(16) 

(17) 

(18) 


Once again these arguments can easily be generalized to an arbitrary initial distribution p{x o). 


3. N immortal or mortal creepers 

Finally, suppose that instead of a single creeper we now have a collection of N statistically independent immortal 
creepers all starting at the same location x$ at time t = 0. The survival probability of the target up to 

time t is simply the probability that it has not been detected by any of the searchers up to that time. We now have 
Sjv(t) = [iS'(t)] JV , and also the asymptotic result SV.oo = <5^ . The same reasoning applies to mortal creepers provided 
they all have the same mortality function, i.e., x = (S^) N . 

III. EXPLICIT RESULTS IN ONE DIMENSION 

In this section we obtain a variety of explicit results for the survival probability of a target in the presence of one or 
more mortal and immortal creepers in dimension d = 1, and compare a number of them with the results of numerical 
simulations. 
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A. Infinite domain 


1. Single immortal creeper 


We begin with a single creeper and a target at location xt g . The convolution structure of Eqs. (JTJ and © translates 
into a Fourier-Laplace transform p(k, s|xo) = J- {p(x, f |xo)} of the propagator which satisfies the so-called Montroll- 
Weiss equation (see e.g. 0,P. 112): 


p(k,s\x 0 ) 


(j)(k,s)e lkx ° 
1 — 4 '(k, s) 


(19) 


We adopt the convention f(k) = J dxf(x)e lkx for the Fourier transform, using the same letter f(k) for the Fourier 
transform as we do for the space-dependent function f(x), again relying on the argument to make the distinction. In 
the particular case of our Markovian walk the above equation takes the form [see Eq. Q] 


p(k,s\x 0 ) 


oo~ 1 e~ ikx ° 

4s) — 1 


( 20 ) 


Next, taking the Fourier-Laplace transform of the joint distribution \Ii(x, i) defined by Eq. (fiTjl. one obtains 


*(k,s) = Z 


+ 


(s + w)l 


2 [ s + w + ikvf s + w — ikvf 
Inserting this expression into Eq. (EU1) yields 


(s + w) 2 + v'jk 2 


( 21 ) 


p(k,s\x 0 ) 


(w + s)e lfcx ° 
s(s + w) + v'jk 2 


( 22 ) 


As is well known, the above expression is simply the Fourier-Laplace transform of the free solution of the telegrapher’s 
equation (see e.g. [27), p. 246): 


d 2 


8 2 


v fQ~2P( x ’ t \ x °) = QpPMvo) 


w dt p ( x,t l x °) 


(23) 


(the telegrapher’s equation can also be obtained directly by differentiating the set of equations H])-© with respect 
to time). Fourier inversion of p(k, s|a;o) yields 


p(x, S\X 0 ) = — J e -^s{s+w)\x-xo\/v f _ ^ 

2vf V s 

Finally, introducing the definition a ffi (s) = y/s(s + w) one has 

p(x, s|x 0 ) = — e -a w ( S )\x-xo\/v f (25) 

2 vf s 

The normalization of this probability density to unity can easily be ascertained through f x Qo p(x, s|xo) dx = s^ 1 . The 
analytic inversion with respect to the Laplace variable s can also be carried out and leads to an expression in terms 
of modified Bessel functions. However, bearing in mind that we are focusing on the survival probability, it is more 
convenient to continue working in Laplace space. 

As an aside, we note that the expression (1241) for the Laplace-transformed propagator can be compared with the 
one for normal diffusion, p n d{x, s|xq). One has 


Pnd{x,s\x 0 ) = £ 


Z ~{x-X 0 ) 2 / [A Dt) 

\ZA-k Dt 


e -yfs/D\x-x o| 

2 VDs 


(26) 


which coincides with the small s limit of Eq. (124[) if one takes D = v'jw 1 . This in turn implies a coincidence at long 
times. 
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The rate at which an immortal creeper detects the target at Xt g any number of times is the sum of the probability 
fluxes of particles that arrive at x tg from the left (with velocity Vf), <& + (x tg ,t), and from the right (with velocity 
—Vf), <&-(xt g ,t). In terms of differential detection probabilities one has 


q(xtg,t\xo)dt = $+{xtg,t)dt + Q-(xt g ,t)dt= / p(x,t-,Vf\xo)dx + / p(x,t;—vf\xo)dx, (27) 




where p(x,t;Vf\xo)dx and p(x,t ; — vj\xo)dx are, respectively, the ensemble probabilities that creepers at a position x 
at time t move toward the right and toward the left. For t > 0 an expansion in dt yields 


q(x t g,t\x 0 ) = Vf \p(x t g,t,Vf\x 0 ) +p{x t g,t,-Vf\x 0 )] = Vfp(x t g,t\x 0 ). 


(28) 


Obviously, q(xt g , 0|xo) = 0 when xo ^ Xt g - However, the case of t = 0 and xq = Xt g must be treated with care. 
First, we realize that p(x, 0; ±v/\xo) = (l/2)S(x — xo) so that, for t —> 0, half of the creepers go to the right from xo, 
whereas none go to xq from the right because there are no creepers beyond xo- In other words, for t —> 0 position Xo 
is special: all the creepers move away from xq, and no creeper goes toward xq. This means that 


g(x 0 , 0 | xq ) dt = $±(x 0 ,0)df = - 


(29) 


However, Vfp{x o, 0|xo)di = <I > + (xo, 0)dt+$_(xo, 0 )dt = 1, so that we can rewrite q(x o, 0|xo) dt = i>/p(xo, 0|xo)di— 1/2. 
On the other hand, p(x tg ,0|xo) = 0 for xo yf x tg . Therefore, these results together with Eq. (1251) can be written as 
follows [28| : 

(Vfp(x tg ,t | X 0 ) , Xo^Xtg 

q(x tg ,t\x 0 ) = i (30) 

[vfp(x tg ,t I x 0 ) - ±5(t), X 0 =x tg 

Note that we have assumed that the walker does not detect the target at t = 0 even if its initial location is Xo = x tg . 
Detection only occurs if the creeper returns to this location. This is not a necessary assumption, and one could easily 
deal with inclusion of immediate target death if the initial creeper location is at the target. 

We close this part by noting that taking the Laplace transform of Eq. (1301) yields a straightforward relation between 
q{xt g , s|xo) and the unrestricted Laplace-transformed propagator p(x, s|xo) of the telegrapher’s equation evaluated at 
x = Xt g . The resulting closed analytic expression can then be used to explicitly compute /(xt s ,s|xo) by means of 
Eq. ©. Inserting the resulting expression into Eq. (flOl) yields an expression for the Laplace-transformed detection 
probability which can subsequently be inverted to explicitly compute S(t). Formally, one may assume that the target 
is killed as soon as it is detected. The detection probability can then be identified with the survival probability of 
the target. If one assumes that the searcher also dies as soon as it detects the target, the target can be regarded as a 
trap and the survival probability of the searcher then also becomes identical to S{t). The computation of the survival 
probability of an immortal walker whose motion is governed by the telegrapher’s equation in the presence of a trap 
was carried out in ref. j29j, and the result is expressible in terms of modified Bessel functions (also see the remarks 
at the end of lHI A3l ). 

2. Single mortal creeper 

We next calculate S^, the asymptotic survival probability of the target in the presence of a mortal creeper. To 
do this we refer back to Eq. (fl8l) . The right hand side requires the Laplace transform of Eq. (l30l) . which is easily 
calculated using Eq. m- We readily obtain 



This can be compared to the corresponding normal diffusion model in the limit Vf —i oo, ui —> oo, with fixed 
D = v'jijj- 1 . In this case one obtains the expression 



l _ e -y/u m /D\x tg -x 0 \ 


(32) 










WWW 

FIG. 1: Plot of the search efficiency in an infinite domain as a function of the mean rate w of direction reversal of the creeper for 
different initial conditions xo and different values of the mortality rate w m : \xt g —xo\ = 0.1 (circles), 1 (squares), 10 (diamonds) 
and 50 (triangles). The solid lines correspond to the analytic result (1331) . whereas the points denoted by symbols are obtained 
numerically from the average over the trajectories of 10 5 independent random searchers. In all cases we use Vf = 1. The 
domain size in the simulations has been chosen large enough for finite-size effects to be negligible. 


for all xo- Note that in the limit of infinite lifetime, w m —> 0, all of these survival probabilities lead (as they should) to 
the result = S* d ^ i Soo = 0, that is, the target does not survive regardless of the initial position of the creeper. 
On the other hand, in the limit of a creeper that dies infinitely quickly, the target survives with unit probability, 
S ^ = S^ d oo = 1, because the creeper essentially does not move before dying. This is the case even when Xq = Xt g 
since our convention is not to count the initial position as a step in the process. 

Comparisons of our results with numerical simulations are presented in terms of the search efficiency 1 — S^, which 
we write explicitly because this makes the comparison more straightforward: 

{ Q^(^to) e —a u (io m )\xtg—xa\/vf 

(tXrri ) + W m 

Gw(w m ) W m 

C^w(w m ) T W m 

Figure |T| shows three plots of the search efficiency as a function of the mean frequency w of change of direction of the 
creeper (direction reversal in dimension d = 1) for different values of the creeper mortality rate w m and for different 
initial conditions. The analytic result (l33l) is in excellent agreement with results from Monte Carlo simulations 
averaged over 10 5 realizations of the random walk. 

As would be expected, the search efficiency increases with decreasing values of the initial distance | Xt g — Xo| and 
decreasing mortality rate (in the latter case, the searcher obviously has a longer time to detect the target). A number 
of the plots in Figure [T] show an optimal value of w for which the search efficiency becomes maximal. This is the 
result of a trade-off between maintaining persistence long enough to avoid spending a great deal of time in a region 
devoid of the target, and breaking persistence often enough to avoid large departures from a nearby target. Obviously, 
the details on how this trade-off works depend on the characteristic time w^ 1 between turns, on the mean lifetime 
w” 1 of the creeper, and on the typical time tq = |xo — Xt g \/vf it would take the creeper to cover the initial distance 
ballistically. In a finite domain of linear size L an additional time scale is introduced, namely the typical time needed 
to cover the full system length in the absence of the target L/vf. This is thus a complex problem, but Figure 1 shows 
the fairly universal feature of an optimal search efficiency (minimum survival probability of the target). 

Figure [2] shows the optimal frequency of turns (optimal persistence) as a function of creeper mortality rate for 
various values of L. While we leave most of the discussion of a finite domain to subsection IIII Bl it is helpful to 
consider the L — > oo behavior shown in Figure [2] For very large mortality rates it is clear that the creeper must 
perform ballistic motion in order to have a chance to detect the target. However, as the mortality rate becomes 
smaller, it may be desirable to change direction a few times in order to avoid long excursions away from the target, 
thereby leading to a finite optimal frequency of turns w op t- The relation between u) opt and w m can then be expressed 
as follows: 

UJ m = (1 ~ W 0 pt|Xtg — Xp\/Vf) 2 

m \x tg - xo\/vf (2 - w op t|x tff - xo\/vf) 

From here we see that w opt = Vf/\x tg — Xo| when w m —> 0, and also that the threshold value of the mortality rate for 
which w opt becomes zero is w m = Vf/(2\x tg — xo|), in agreement with Figure [2] and with the description given above. 
We return to Figure [2] later to discuss the case of finite L. 


x 0 Xt g 

Xq — Xfg 


(33) 
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FIG. 2: Behavior of the optimal persistence w op t as a function of the mortality rate u i m for different values of system size L, 
L = 100 (circles), 10 3 (squares), 10 4 (diamonds) and L —> oo (triangles). Solid lines correspond to the analytic result (1601) . 
while symbols denote data points obtained from Monte Carlo simulations (averaging over the trajectories of 10 5 independent 
mortal creepers). In all cases we place the target at xt g = 0, begin the walk of the creeper at xo = 10 and the creeper has 
speed Vf = 1. 


3. Time dependence 

Results for the search efficiency 1 — S*(t) at finite times can be obtained by inversion of the exact expression in the 
Laplace domain. The latter reads 


1 


(s + UJ + LUm ) 1 / 2 


1 . J s (s + Wm ) 1 / 2 + (s + UJ + Wm ) 1 / 2 

S I 1 (s + UJ + Urn) 1 ! 2 ~ (s + UJm) 1 ! 2 

S (s + UJ + Urn ) 1 / 2 + (s + UJ m ) 1 / 2 ' 


o-<x u (s+u m )\x tg -x 0 \ /v /) Xq ^ Xtg 


(35) 


Xq — 


Let us first focus on the case where the searcher starts at the target location (xo = Xt g ). Inverting the equation for 
the Laplace transform by means of the Faltung theorem and the shift theorem, one obtains 


1 — S*[t) = J* dr e~ UmT C~\ T ^ 


((S + w ) 1 / 2 — S 1 / 2 ) 2 


UJ 


According to ref. [3(j (p. 215, formula 79), one has 


C~ l = 


((s + a;) 1 / 2 — s 1 / 2 ) 2 




where Ii(-) stands for a modified Bessel function. We thus find 


(36) 


(37) 


1 - S*(t) = dTT- 1 e~^ +2uim)T/2 I i (y) . 


(38) 


One can easily check with look-up tables that the value of the integral on the right hand side in the limit t —> oo 
coincides with the value of 1 — S ^ given in Eq. (l33l) . In passing, we also note that in the absence of mortality 
(w m = 0) the integral is known for all times t (see e.g. ref. [3l|, p. 46). One finds 


! - S(t) = 1 - e~^ 2 h (y) - e-^Io (y ) . 


( 39 ) 


Returning to the general case u> m ^ 0, the killing probability can be rewritten as follows: 
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1 - S*(t) = 


&OJ i^m ) ^77 

Q-oj (^m) H - Wn 


j™ dr T -'e-^ +2 ^ T ' 2 h (p . 


For long (but finite) times t , we now use in the above formula the asymptotic expansion for large x: 

{l-± 15 


h(x) = 


o 0 


-3> 


p2irx l Sx 128x 2 

Neglecting all the terms beyond the first subdominant term and making use of partial integration we find 


J°° dTT-^+^^h (p 


( UT \, 

1 

/»oo 

V 2 ) ' 

~ v^J 

( 


dr: 


(wr) 1 / 2 4^/n 

3 


dr ' 


t 

—OJm t 


(lot) 3 / 2 

— UJmt 

+ o'" 


TTUJ U) m t 3 / 2 \frLO (2W m + 4w)w m t 5 / 2 \ t 7 / 2 


At long times one thus has 


1 ~S*(t) = 


LT uj (oJ rr i ) LO n 


1 e~ Umt 3 

+ 


(co m ) + W m i/™ oo m t 3 / 2 a/7tw (2w m + 4w)cc m t 5 / 2 
In the opposite limit of early times we can make use of the Taylor expansion 


O 


tV 2 


(40) 


(41) 


(42) 


(43) 


h( x ) = P^ + 0(x 5 ) 

and the series definition of the exponential function in (1381) to obtain 

co T 2cu 


1 -»•<*)-? i dT[1 -(^) T+ ( 


r+( g + (u i± 2a^') r2 + 0 ( T 3, 


Thus, 


1 w(w + 2w m ) i2 ( f co 3 , w(w + 2w m ) 2 ^ 3 , 

1-6 (t)- T jg t + \J84 + 96 J + 


(44) 


(45) 


(46) 


We now turn to the case x$ ^ x tg . Taking the inverse Laplace transform of the expression given in (1351) for this 
case, one finds 


1 - S*(t) = to- 1 J* dr e-^ r C~X T {(a + cn - a u (s))e - a “^ T °} , (47) 

where we have introduced tq = vj 1 1 x tg — a?o | ■ In what follows we restrict ourselves to the case t > r 0 , since for t < r 0 
one trivially has 1 — S*(t) = 0. Now we invoke the following results given in ref. [32| : 


CjXr [ e ~ a ^ 

C-\ T {se~ a ^ s)T0 

^ {«<-(«) e- a “ (s)7n 
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where 8(-) stands for the Heaviside step function. Substituting these results into Eq. (l47l) . applying partial integration 
for integrals involving derivatives of Bessel functions and carrying out integrals not involving Bessel functions, we are 
left with the following expression: 



where additional use of the identity ( d/dx)Io{x ) = I\(x) has been made to express the time derivative of the zeroth- 
order modified Bessel function in terms of the first order modified Bessel function. 

Note that there is a discontinuity of 1 — S*[t) at t = to, where it jumps from zero to a finite value. In this case, 
1 — S'* (to) is the probability that the target is found by the creeper exactly at t = tq. For this to happen, i ) the 
creeper must start moving towards the target at time to, ii) it must not change direction in the time interval between 
0 and To and, in) it must not not die during this time interval. The probabilities for these three independent events 
are respectively 1/2, e - “ T °/ 2 and e~ ulmT ° . The product of the three yields the result 1 — S*(t 0 ) = (1/2) e-K 2 ”")^/ 2 , 

Let us now focus on the study of the asymptotic behavior of (l5lT) . We first start with the long time regime. Making 
use of the x-large expansion 


Io(x) 
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128 a ; 2 


+ 0(x 



(52) 


together with (TUI) , and applying partial integration to the corresponding integrals, one finds the long time behavior: 


x _ s *, f j = _(LW_ e -CL^ m )T 0 _ { l + (u;+^m)To \ e 

a u (u m )+w ra V 2u] m V™ ) t 3 / 2 

f3(2 uj + LUm)(l + (w + w m )r 0 )\ e -aJmt 

V 8vW< J V2 + 1, t V2 J ' 

Let us now consider the opposite limit of short times. Here, one uses 


( 53 ) 


I 0 ( x ) = l + ^ + ^+O(x e ) (54) 

together with (14411 to perform an expansion of the functions appearing in (1511) in the small parameter A = t — tq. The 
final result is 


1 - S*(t) = e -(“+ 2 “m)ro/ 2 + w(wto_+ 2 ) A + (euro - 4 )ui 3 tq - 16 UJ 2 - 16 (o;to + 2 )uuj m + Q ^ A 3 ^j 


(55) 


Note that the linear term in A does not depend on ui m . The mortality rate only appears in the quadratic order 
term. 

Finally, we note that in the absence of mortality, the formula (l5lT) takes the simplified form 


1 -S( t ) =*-*”»(| 



( 56 ) 


Alternatively, this formula can be obtained as the spatial integral of the propagator given in ref. [29j for the telegra¬ 
pher’s equation in the presence of a single trap. 
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FIG. 3: Log-linear plot of the cumulative probability to find the target as a function of time. The left figure corresponds to 
the case Xt g = xo, whereas the right figure corresponds to the case | xt g — xo\ = 10. In both cases we take ca m = 0.002. The 
symbols correspond to simulation results for to = 0.01 (circles), to = 0.1 (squares), to = 1 (diamonds) and to = 10 (triangles). 
The solid lines correspond to the analytic results respectively given by the formulae (1381) and «5D- 


Returning to the case with mortality, Fig. [3] shows a comparison for different parameter values between analytical 
and numerical results for the detection probability as a function of time. Excellent agreement is found. 

We close by noting that, in the presence of creeper mortality, the survival probability of the target S*(t) and that 
of the creeper S* r (t) become different if one assumes that both creeper and target die instantaneously upon detection. 
However, there is a straightforward relation between both of them, namely, S*(t) = S* r (t) + oo m fg S* r (r)dT (see ref. 
in for a detailed derivation). Obviously, one must have S*(t) > S* r (t) in this case, since the creeper can die not only 
upon encounter with the target, but also spontaneously. 


It- N mortal creepers 

Consider the case of TV > 1 independent creepers that begin their walk at the same location and at the same time. 
The question in this case is whether an optimal search efficiency of the walkers (i.e., minimum value of the survival 
probability of the target) is still found. The search efficiency is now simply 1 — Sn,oo = 1 — (S* a ) N . The condition 
for w op t is again found by requiring that the derivative of the efficiency with respect to to be zero. However, the 
condition for this derivative to become zero is the same as for the TV = 1 case. Therefore, the value of w op t is the 
same. However, since, S ^ < 1, the detection probability of the target tends to one as TV —> oo. In words, when the 
number of walkers becomes infinitely large, one of them will sooner or later reach the target with certainty no matter 
how quickly the walkers die on average. 


B. A finite domain 

1. Single mortal creeper 

The most straightforward extension of our results to a finite system is to a ring of length L. On this ring, we place 
a point target and calculate its survival probability S^ in the presence of a single mortal creeper. This is most easily 
done by exploiting the relation between the probability density p{x,t\xa) introduced in Eq. © and that of the finite 
system p L (x,t\x 0 ): 


OO 

p L {x,t\x 0 ) = E p(x + nL, t\xo), 0 < x < L. 

n =—oo 


In Laplace space, this translates into the relation 

OO 

p L (x,s\x 0 ) = E p(x + nL , s|rco), 0 < x < L. 


(57) 


n =—oo 


(58) 
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FIG. 4: Plot of the search efficiency in a finite domain (L = 200) as a function of uj for different values of the mortality rate and 
for different initial conditions: \xt g — xo| = 0.1 (circles), 1 (squares), 10 (diamonds), and 50 (triangles). Solid lines correspond 
to the analytic result (16011 and symbols are obtained numerically from the average over the trajectories of 10 5 independent 
random creepers. In all cases we take Vf = 1. 


Using Eq. (1251) we then obtain an expression in terms of a geometric series which can be summed. The sum yields 


p L (x,s\x 0 ) 


1 tru;(s) 
2 vf s 


e -a w (s) \x-x 0 \/v f _|_ e -a m (s) (L-\x-x 0 \)/v s 
1 _ g -a w (s) L/vf 


0 < x < L. 


(59) 


One can easily check that the above equation is normalized, i.e., dxpL^x, s |0) = s -1 . The limit of an infinite 

system is also recovered as one lets L —> oo in the above expression. 

The corresponding detection rate qL(xt g ,t\xo) is obtained by making use of Eq. (|59l) in Eq. (l30l) . The survival 
probability of the target is subsequently computed from Eq. (El), which remains valid for a system of finite size. The 
final result reads as follows: 


1 - 


j (u-'rn ) 


Q* _ 

^OO 


^rnP[ T O:^; (cU m )/3^ 


( e -au(um)\xtg-Xo\/vf _j_ e ~ a “ (ui m ) (L-\X t g ~X 0 \)/vf ) ^ X(J ^ 


P [ 


(60) 


. ^mPj^ T 


•t-0 — tgi 


With Pf = 1 ± e -a^rn.)L/v f . 

We note the great qualitative and even quantitative similarity between the results for L finite (Figure 0]) and those 
in Figure [l] This similarity is also seen quite explicitly when comparing Eqs. (1601) and (1311) in the limits aj m —> 0 and 
com —> oo. The differences between the figures and equations of course arise from our domain size (L = 200) which 
is small enough to produce non-negligible finite size effects in the ranges of values we chose for the other parameters. 
Without loss of generality we place the target at x tg = 0 or, equivalently, x tg = L. A physical sense of the situation 
is obtained as follows. Suppose that the creeper begins its walk at a location near the target (the greatest initial 
distance between creeper and target occurs when the creeper begins its motion near xq = L/ 2). If the creeper moves 
toward the target along the shorter direction without a direction reversal, it reaches the target very quickly and the 
target does not survive these events, thus decreasing the survival probability S^. These events are independent of 
system size, and thus are not affected by the finite size of our system. If the creeper moves toward the target along 
the long side, the minimum time needed for the creeper to reach the target at xt g = L is of order L/vf. The mean 
time to death of the creeper is of order l/w m . If this time is shorter than L/vf , the creeper dies before reaching the 
target, and finite system size effects are again not seen on average. This description fits that shown in the right panel 
of Figure U since L/vf = 200 for the parameter choices indicated in the caption, and 1 /uj m = 50. If the creeper 
suffers many trajectory reversals, then it is even more likely that it dies before reaching the target. The left panel is 
the one in which the parameter choices would most clearly lead to finite size effects since now l/ui m = 5000 and for 
sufficiently low reversal frequency the searcher will almost certainly reach the target, even though its initial position 
is not so close to the target. The survival probability of the target S ^ should be smaller than in the right panel, and 
this is what we see. However, the results for the infinite system in Figure [Hand for the finite system in Figure H] are 
rather similar, not only qualitatively but even quantitatively for all the parameters considered. 

Earlier we introduced Figured which shows the reversal frequency uj of the creeper trajectory at which the search 
efficiency 1 — S ^ is a maximum as a function of the mortality rate of the creeper. We called this frequency the optimal 
persistence w = u; op t. We see there that the optimal persistence depends quite sensitively on L in some regimes. For 
very small systems the optimal value is always near zero (as happens for L = 100 in Figure [2]), that is, regardless of 
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FIG. 5: Schematic representation of the three regimes observed in Figure [2] for intermediate values of L. Sketches of typical 
trajectories showing how the creeper (the black point) most likely reaches the target are shown as arrows for each region. The 
target is located at Xt g = 0, which is equivalent to Xt g = L because of the periodic boundary conditions. 


the mortality rate of the walker the best strategy is ballistic motion. As L increases, a new regime emerges where 
w op t is non-zero, that is, the creeper is most likely to step on the target with reversals in the direction of motion. 
For very high creeper mortality rate, the best strategy is again ballistic motion. This is consistent with the behavior 
seen in Figure 2) Note that the maximum value of w opt is roughly the same for all the L’s shown in Figure [2] but 
the width of w m -interval where w 0 pt is appreciably different from zero is strongly L-dependent and it is narrower for 
smaller L. Finally, as L —> oo the region with a finite optimal persistence increases in the direction of low mortality 
and eventually becomes dominant. 

In order to further clarify these results, Figure [5] shows a sketch of typical trajectories of the creeper in the three 
regimes observed for intermediate values of L. For a high mortality rate (Region III), the only way the searcher 
can successfully hit the target is to move ballistically toward it along the shortest path before dying, and therefore 
maximum persistence is desirable. However, as the mortality rate decreases, the creeper has time to reverse its 
trajectory and move in one direction or the other, eventually with equal probability. Introducing some breaking of 
persistence in this regime is seen to be a desirable strategy for the creeper. The reason is that in those runs where 
the creeper starts moving away from the target, a change in direction makes it more likely to detect the target before 
dying. Finally, if mortality is so low (Region I) that the searcher has enough time to reach the target location x = L 
by moving ballistically before dying, one again finds that uj op t = 0. Breaking the persistence in this regime would 
evidently introduce unnecessary overlapping of the trajectory with itself and increase the cumulative probability of 
death due to the additional delay. As one would expect, in the limit L —> oo Region II grows at the expense of region 
I (which becomes vanishingly small), and uiopt saturates at a value that depends on the initial separation |a;o| (recall 
that we have set Xt g = 0 for simplicity). 

From our analysis so far, it is clear that the transition from Region II to Region III points to the fact that the 
characteristic lifetime u;” 1 of the creeper has gone from being longer to being shorter than xo /v/, while the transition 
from Region II to Region I occurs as w” 1 grows and becomes comparable to L/vf , the typical time it takes the creeper 
to reach the target if the motion is ballistic. Implicit in this reasoning is the condition Xq <C L/ 2; otherwise, the 
necessary trade-off for the onset of an optimum persistence breaks down. This is the reason why in Figure [2] no onset 
of Region II is observed for sufficiently small values of L. 

Mathematically, the specific condition that ensures the existence of Region II so that the search efficiency is not 
always maximal for w = 0 arises from the requirement c?(l — S 00 *)/dui\ ul -o > 0. Using Eq. (16U1) in this requirement 
one finds that the condition 

(l-x o + x 0 e tJm(i “ 2a:o/ ^ ) ) + e ~ 2urnX °lv/ + e ~ u ™ L / v f - e ^(L-2x 0 )/v f _ 1 < Q 


Vf 


(61) 
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co 

FIG. 6: Asymptotic span of a Id mortal creeper as a function of the frequency of reorientation for different mortality rates (we 
have set Vf = 1). 


must be fulfilled. From this condition one can establish that both uj m L/vf and Xo -C L/2 must necessarily hold 
for Region II to exist. Summarizing, a non-zero optimal persistence can only arise if the searcher is initially within a 
short distance of the target. Distant targets can not be reached within the creeper’s lifetime. 


2. N mortal creepers and the target problem 


Next we again consider a target in the presence of N creepers. If the creepers start from the same location, the 
remarks made in the case of an infinite system apply here as well. More interesting is the situation where the position 
of each creeper is drawn from a uniform distribution p[x o) = L~ l in the interval [—L/2, L/2], with the target located 
at the center of the interval (xt g = 0). This so-called target problem has long been discussed in the literature for 
immortal random walkers and more recently by three of us for mortal diffusing particles and mortal walkers with 
independent stepping time and stepping length distributions EM. The survival probability of the target at 
time t in the presence of the mortal creepers is given by the same expression as for usual walkers (leapers), namely, 


/ ^ j-L /2 \ / 2 r L / 2 

S%{t)= l-j dx 0 S*(t\x 0 )J = ( -j dx 0 S*{t\x 0 ) 


(62) 


In the thermodynamic limit we let both L and N go to infinity, keeping the initial density of mortal creepers fixed at 
a value p o = N/L. The above equation then becomes 

2p 0 r N/{2po) 

~N~ Jo 

where S*(t) denotes the survival probability of the target in the aforementioned thermodynamic limit. For t — > oo, 
we now use Eq. (1531) to evaluate the last integral explicitly: 

Sla = lim S*{t) = exp{—2 p 0 Vf/{u} m + a w (w m ))}. (64) 


S*(t) = lim Str(t) = lim 

N-too N-> oo 


dx o (1 


N 


= exp{-2p 0 / dx 0 (1 - S* (£))}, (63) 


Note that the corresponding search efficiency 1 — decreases steadily with increasing average direction reversal 
frequency w, implying that no optimal persistence w op t > 0 arises in this limit. 

There is an interesting and useful connection between the last term of Eq. (|63l) and the territory explored by a 
single walker during time t Explored territory includes all those locations visited any number of times. 

For mortal creepers, this territory is given by the integral in the argument of the exponential in Eq. (1631) . that is, 
2 J 0 °° dx o (1 — S*(t)). Thus, according to Eq. (l64l) . in the limit t —> oo the territory explored by a single mortal creeper 
is 2 Vf/(u> m + a u (w m )) (in Id this quantity is identical with the span of the mortal random walk in the absence of the 
target). This analytic result is confirmed by numerical simulations (see Fig. [6]). 





16 



FIG. 7: Optimal persistence obtained as a function of the mortality rate cu m for periodic square lattices of size LxL, with L = 25 
(circles), 100 (squares), 400 (diamonds) and L — > oo (triangles). The results were obtained by averaging over trajectories of 10 s 
independent mortal creepers. The target is assumed to be a circle with unit radius centered at the origin (xt = ( xt g , yt) = 0). 
In all cases the parameter values \J ( xtg — xo) 2 + (yt — yo) 2 = 2 and Vf = 1 were taken. 


IV. RESULTS IN TWO DIMENSIONS 

In order to assess whether the onset of an optimal persistence is specific to the one-dimensional case and the pecu¬ 
liarities of space exploration properties of random creepers in this constrained geometry, we have briefly investigated 
the behavior of the search efficiency in a two-dimensional system. The route to an analytic solution appears to be 
cumbersome, but numerical simulations allow us to conclude that an optimal persistence also exists in this case. 

We again assume perfect detection, implying that the target will instantaneously be detected whenever the trajectory 
of a pointwise searcher attempts to penetrate it. Turns interrupting the ballistic trajectories of the searchers occur 
according to a uniform turn angle distribution. In particular, this also applies for the velocity initial condition. 

Results for the search efficiency in two-dimensional torii (arising from periodic boundary conditions) are displayed 
in Figure [3 Trends similar to those of the one-dimensional case displayed in Figure [2] are observed. Once again, 
regions I, II and III can be clearly identified and the behavior of the optimal persistence as a function of L follows the 
same trends as in dimension d = 1. A reasoning similar to the one underlying Figures [2] and [5] is expected to apply 
here as well. 


V. SUMMARY AND OUTLOOK 

We have investigated the effect of mortality processes in the efficiency of a search process performed by a class of 
randomly moving particles with spatiotemporal coupling. Our main conclusion is that for a suitable parameter choice 
a finite optimal persistence maximizing the long-time detection probability of the target appears. This optimization 
effect has been shown to exist both in finite and infinite Id domains, and to persist in the 2d case (in the case of an 
infinite one-dimensional domain we have also derived an explicit expression for the survival probability of the target 
after a finite time t). However, optimization with respect to the frequency of turns is suppressed in the limit where 
the target is surrounded by a finite density of mortal creepers initially distributed at random over an infinite interval 
(Id target problem). 

Even though our model is a strong simplification of real systems (say, a predator hunting a prey), the robustness 
of the aforementioned effect leads us to believe that it may also play a major role in more realistic systems. In any 
case, our model incorporates features which are likely to be relevant in an ecological context, namely, the fact that 
long jumps must pay a time penalty, and lifetimes which are comparable to the relevant time scales of transport [28| . 

The present work can be extended in many interesting ways. Further efforts are indeed needed to completely 
ascertain the role of spatiotemporal coupling and mortality in a broad class of target search problems. As already 
mentioned in the introduction, one could for instance replace the exponentially decreasing waiting time density with 
a long-tailed waiting time density and investigate whether the resulting walk leads to similar optimization effects as 
the ones observed here. Such problems are interesting even in the absence of mortality. Bearing in mind possible 
applications in ecology, the case of Levy flights with long-tailed jump length distributions but no coupling between 
jump lengths and waiting times could be investigated (in this case special care should be taken when defining the 
details of the detection process). Finally, one could relax the condition of an immobile target and assess the behavior 
of the search efficiency by devising suitable approximations for the case where the target is also allowed to move. The 
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question is whether some kind of Pascal principle applies, as in the case of conventional Brownian motion 
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